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The criteria for the existence of supersonic and multi- 
ple topological excitations (kinks) in the driven Frenkel- 
Kontorova model (a chain of atoms placed into an external pe- 
riodic potential) with anharmonic (exponential) interatomic 
interactions are studied. 



I. INTRODUCTION 

Topological excitations play a very important role 
in system dynamics, because they are responsible for 
mass and charge transport in solids and on crystal sur- 
faces. As classical examples, one can mention disloca- 
tions and crowdions described by the Frenkel-Kontorova 
(FK) model yj, where the topological excitations cor- 
respond to kinks that describe local compression or ex- 
pansion of a commensurate structure. The FK model has 
numerous applications in superionic conductivity, surface 
physics, hydrogen-bonded chains, Josephson junctions, 
tribology, etc. (e.g., see |2| and references therein). 

In the continuum limit approximation, the equation of 
motion of the FK model reduces to the exactly integrable 
sine-Gordon (SG) equation. But in continuum models, 
even in a model with anharmonic (but local) interac- 
tion, the topological excitations are always subsonic, the 
kink cannot propagate with a velocity v larger than the 
sound speed c because of Lorentz contraction of kink's 
width. Moreover, in the classical FK model, the kinks 
of the same topological charge repel from one another 
and, therefore, they cannot carry a multiple topological 
charge. 

However, simulation demonstrates that supersonic 
kinks as well as multiple kinks do exist. For example, 
Fig. [j] shows the propagation of supersonic single and 
double kinks in the FK model with exponential inter- 
atomic interaction. The motion of topological solitons 
with supersonic velocities was firstly predicted, to the 
best of our knowledge, analytically by Kosevich and Ko- 
valev H in the FK model with some specific interatomic 
interaction in the continuum limit approximation. Later, 
supersonic topological solitons were observed by Bishop 
et al. Q] in molecular dynamics study of polyacetylene. 
Then the supersonic kinks were studied numerically in 
the discrete FK model with anharmonic interaction by 
Savin §. It was shown that for certain supersonic kink 



velocities, when its width coincides with that of the cor- 
responding Toda soliton [0], the kink propagates almost 
without energy losses. Multiple fast (but subsonic) kinks 
were firstly observed numerically by Peyrard and Kruskal 
[0 in the classical highly discrete FK model. Alfimov et 
al. |^| have shown that multiple kinks exist also in con- 
tinuum systems with nonlocal interaction. The bounded 
states of kinks (subsonic as well as supersonic) in the case 
of anharmonic interaction were also studied numerically 
by Zolotaryuk et al. ||. It was found that these multi- 
ple kinks are asymptotically unstable. The dynamics of 
the generalized FK chain driven by a dc external force 
was studied numerically in [ p~0|JTT| , where we observed the 
existence of supersonic kinks and multiple (double and 
triple at least) kinks. Recently, the existence of multiple 
kinks for certain kink velocities in the discrete FK-type 
model was proven rigorously . 
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FIG. 1. Atomic trajectories of the FK model with expo- 
nential interaction for /3 = l/V, g = 1, and n — 0.05. (a) The 
single supersonic kink, / = 0.45, v^/c w 1.28, and (b) the 
double kink, / = 0.60, v^/c « 1.75. 
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The aim of the present paper is to find the criteria for 
the existence of supersonic and multiple topological exci- 
tations in the FK-type models. We will show that, first, 
the model must be discrete as was already mentioned 
above. Second, because kink's motion in a discrete chain 
is damped due to radiation of phonons, we must apply 
a driving force to support the motion (this point was 
lost in the previous studies |§-[|,0U). Third, the in- 
teratomic interaction must be anharmonic. Under these 
three conditions, the model admits both supersonic kinks 
and multiple kinks. 

The paper is organized as follows. The model is intro- 
duced in Sec. ||. Then in Sec. Ill the problem is studied 
with the help of the continuum limit approximation. We 
show that if the discreteness effects are properly taken 
into account, the model formally allows both supersonic 
kinks and multiple kinks. Moreover, we prove that for 
a given set of parameters the model admits either the 
single-kink solution or the double-kink solution, but not 
both solutions simultaneously. Numerical solution of the 
corresponding ordinary differential equation, which can 
be obtained with any desired accuracy, does show the 



existence of supersonic and multiple kinks. In Sees. IV 
and ^ we develop an approximate variational approach 
that helps to find the conditions for existence of super- 
sonic and multiple kinks. We show that supersonic kinks 
may exist in the model with anharmonic interaction only, 
and the multiple kinks may be stable for supersonic ve- 
locities only. These conclusions are confirmed by sim- 
ulation. Finally, Sec. VI concludes the paper with the 



comment that supersonic kinks and multiple kinks may 
be considered as "disturbed" Toda solitons. 



II. MODEL 

We consider the chain of atoms placed into the sinu- 
soidal substrate potential and driven by a dc force / ap- 
plied to all atoms, so that the equation of motion is 
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xi{t) + nxi(t) + — [V(xi+i - xi) + V(xi - xi-i) 
+ sin xi(t) - f = 0, 



(1) 



where xi is the coordinate of the Zth atom, r\ is the exter- 
nal viscous damping coefficient introduced to compensate 
the driving force, and 



V(x) = V Q exp(-0x) 



(2) 



describes the exponential repulsion between nearest 
neighboring atoms. We assume periodic boundary con- 
ditions with the number of atoms N — N s + p, where 
N s is the number of wells of the periodic potential, so 
we have one multiple p-kink (jp excessive atoms) inserted 
into the commensurate structure. Similarly to the Toda 



model ||, one may add also an attractive linear inter- 
action, so that the interaction reduces to the harmonic 
one in the limit (3 — ► (the classical FK model) and to 
the hard-core interaction in the opposite limit (3 — > oo. 
Throughout the paper, we use the dimensionless system 
of units, where the atomic mass is m = 1, the period 
of the external potential is a = 2tt and its amplitude 
is e = 2. In these units, a characteristic frequency of 
atomic vibration at a minimum of the external potential 
is luq = 1, a characteristic time scale is tq = 2ir, and 
the maximum value of the external dc force, when the 
minima of the sinusoidal substrate potential disappear 
and the topological excitations cannot exist anymore, is 
/=!• 



III. CONTINUUM APPROXIMATION 

The systematic procedure to derive the equation of mo- 
tion in the continuum limit starting from the discrete lat- 
tice was proposed by Rosenau |0|| . For the anharmonic 
FK model (|]) it leads, in the first order of the discreteness 
parameter, to the equation 



u+r]u — d u" (1 — adu') + sin u — f 
-h 2 u" + (it') 2 sin it — u" cosit = 0. 



(3) 



where d = a^/g is the width of the static kink, g is the 
elastic constant defined as g = V"(a) — Vof3 2 exp(— (3a), 
the anharmonicity parameter a is defined as 



a V'"(a) _ J_ 
dV"(a) ~ y/g' 



(4) 



and the parameter h 2 — a 2 / '12 = 7r 2 /3 describes the dis- 
creteness effects. Looking for a traveling- wave solution 
of Eq. (|^) in the form u(x, t) — u^{x — vt) = u^(z), we 
obtain the equation 

h 2 v 2 u'k" + (c 2 - v 2 - h 2 cosu k ) tt£ + h 2 (u' k ) 2 sinuk 

— aePuk'itk + rjvu'^ — sinitk + / = 0, (5) 

where c = 2n^/g is the sound speed (in our system of 
units c = d). 

Although the travelling-wave ansatz may be too crude 
due to radiation of phonons by the moving kink (e.g., see 
Fig. |l|a), Eq. (||) allows to find kink's asymptotic rigor- 
ously, because the radiation has to decay due to nonzero 
damping coefficient r\ in the model under consideration. 
At z — > oo, substituting u^(z) — Uf oc exp(— z/d±) into 
Eq. (H), we obtain for the kink width di(v) the equation 



d\ cos Uf + d\nv — d\ (c 2 



h 2 cosu f ) = (/iw) 2 , (6) 



1/2 

where cos it/ = (1 — / 2 ) . Similarly, we can find the 
tail asymptotic behind the kink, itk(z) — (it/ + 2irp) oc 
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exp(z/c?2) at z — > — oo; the width d 2 (y) has to satisfy the 
same equation (||) but with v — ► — u. One can see that at 
low velocities, \v\ <c c, the discreteness effects lead to a 
decrease of the kink width in agreement with theory Q . 
However, now Eq. (^|) has a solution for any kink velocity 
v . Thus, the discreteness effects remove the restriction 
\v\ < c of SG-type equations. Now even for the classical 
FK model with harmonic interaction the kink may move 
with any velocity v. 

A kink solution corresponds to a separatrix of the con- 
tinuum limit equation. To find the separatrix of Eq. (|^), 
let us normalize the coordinate z = z/d, the velocity 
v = v/c, and define the dimensionless discreteness pa- 
rameter h = h/d = l/^/12g. Introducing the new vari- 
able £ = Uk(?) and the function w(t;) = u' k (z), Eq. (||) 
can be rewritten as 



j/iV [«/"(0™ 2 (0 + W(£K(£M£) + K(£)f 

- aw'(£)w(Z) + (l - v 2 - h 2 cos^j w'(0 

+ (h 2 sin*;) w{£) + - sin£ + / = 0. (7) 

A (multiple) kink solution has to satisfy the bound- 
ary conditions u^{— oo) = Uf + 2np, u^(+co) = Uf, 
u' k (±oo) = 0, or 



w(uf + 2irp) — w(uf) = 0. 



(8) 



For example, if w(£) is a separatrix for the double kink 
(p = 2), it has to connect the points (£ = Uf + 4n, w = 0) 
and (£ = Uf,w = 0). However, due to periodicity of 
the substrate potential the function w(£ + 2tt) must cor- 
respond to the separatrix solution as well. Thus, on the 
(£, w) plane the separatrices of the multiple (p > 2) kinks 
must intersect at some point with w 7^ 0. One can show 
that in the model without the discreteness effects, h = 0, 
when the phase space of E q. (|7| ) is two-dimensional, such 
intersections are forbidden EW . Thus, Eq. (^) with h = 
allows neither supersonic kinks nor multiple kinks. On 
the contrary, at h ^ the phase space of Eq. ([7]) is 
four-dimensional, thus the separatrices corresponded to 
multiple kinks may not intersect, and multiple kinks are 
allowed in principle |l2| . Thus, the discrete model for- 
mally allows both supersonic and multiple topological ex- 
citations. 

Although we cannot find the separatrix solution an- 
alytically, the ordinary differential equation (Q) can be 
solved numerically with any desired accuracy. Thus, if 
one could find a separatrix solution corresponded to su- 
personic or multiple kink, this will prove their existence. 
Indeed, looking for a separatrix solution numerically for 
the P — 1/tt and 77 = 0.05 case, we found that at small 
discreteness, g = 10, so that h w 0.09 and a w 0.1, 
the separatrix solution corresponds to the 2-7r-kink at 
forces as large as / = 0.9, when the kink is supersonic, 



«k/c w 1.13. On the other hand, for higher discrete- 
ness, g — 1 so that h w 0.29, we saw the 27r-kink at 
/ < 0.2 when v k /c < 1, and the 47r-kink at / > 0.6 
when «k/c > 1.3. Thus, both supersonic kinks and multi- 
ple kinks do exist, at least for some particular choices of 
model parameters. 

Near kink's tails, z — > ±00, e.g. for £ = Uf + 2-kji + e, 
where |e| -C 1, we can use the expansion 
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«;(£) = axe + -a 2 e* + -a 3 e° 



(9) 



where a\ = w'(v,f), a 2 = w"(uf), etc. Substituting this 
expansion into Eq. (]7j) and grouping together the terms 
of the same power of e, we obtain the following equation 
for ai [cf. with Eq. (§)], 

T2—2 4 1 /-, ~2 T2 \ 2 1 ~ n 

ft o Bj + (1 - v — h cosu/)a 1 + r/vai — cosu/ = 0, 

(10) 

which always has two roots, one positive and one nega- 
tive, for any kink velocity v. Then, equating the terms for 
higher powers of e, we obtain the relations that uniquely 
determine the coefficients a 2 , 03, etc. Thus, the separa- 
trix solution of Eq. (Q) is unique, i.e., for a given set of 
system parameters the model has either the single-kink 
solution or the double-kink one, but never both solutions 
simultaneously. 

Thus, we have demonstrated the existence of super- 
sonic and multiple kinks in the driven discrete FK model 
with anharmonic interaction (note that the discreteness 
of the model and the anharmonicity of the interaction 
are the necessary conditions). But to find the parame- 
ter range for their existence (the sufficient conditions), 
we have to study numerically either the continuum-limit 
equation (]?]) or, better, the original discrete model (Q). 
The approximate variational approach described below, 
essentially simplifies this task. 



IV. SUPERSONIC KINKS 

It is easy to show that in the case of / = r\ = h = 0, Eq. 
(||) corresponds to an extremum of the following energy 
functional [(TBI, 



E[u{z)\ = dz 



c 2 — v 2 , ,,2 ad 



(u) (u) —cosu 



Substituting a simple SG-type ansatz 

usg{z) = 4tan _1 exp(-z/d c ff) 
into Eq. (|ll]), we obtain 

- 2 - v 2 2tt c 3 



E(d eS ) = 4- 



a +4d eff . 

err 



(11) 

(12) 
(13) 
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Using Eq. (|l6|), Eq. (|l4|) can be rewritten in the form 
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FIG. 2. The velocity Wk of the single kink versus the force 
for ft = 1/tt, g = 1, and r\ = 0.05. The solid curve is for 
simulation results, and the dashed curve, for variational ap- 
proximation. Inset: the critical kink velocity v c (/3) at the 
fixed force / = 0.5. 

Although the variational approach does not describe 
rigorously the kink tails because of neglecting the dis- 
creteness effects, it allows us to find analytically the 
shape of the kink's core and, therefore, to calculate ap- 
proximately the kink velocity for the model with anhar- 
monic interaction. Indeed, looking for extrema of the 
function E(d c g), we come to the equation E'{d e g) = 0, 
or 



« 3 = 



1 



7T 

k + -a, 



(14) 



where we introduced the new variable n — d e g/d. For the 
harmonic interaction, a = 0, Eq. ( |l4| ) has a solution for 
\v\ < c only, which describes relativistic narrowing of the 

1/2 

SG kink, k = [l — (v/c) 2 ] . But for the anharmonic 
interaction, a > 0, Eq. ( |l4| ) has a solution for any kink 
velocity v, including supersonic velocities \v\ > c. We em- 
phasize that supersonic excitations are possible for kinks 
(local compressions) only. 

Considering the kink as a rigid quasiparticle, the kink 
effective mass can be introduced as (e.g., see j§) 



1 



m k = 



dz [u'(z)Y 



nd e a 



(15) 



Then, assuming that kink's parameters at nonzero / and 
rj are the same as those for the / = r\ = case, the 
steady-state kink velocity can be found approximately 
from the equation 



Vk = f/nikV = TrcfK(v k )/4i]. 



(16) 



, *7 
i+ I -r 

47/ 



7T 

k + -a. 



(17) 



Numerical solution of Eq. (|l7j ) allows us to find the func- 
tion w k val '(/) that is shown by the dashed curve in Fig. 
H together with the dependence v^if) obtained by solu- 
tion of the discrete equation of motion (fy) . One can see 



(var) 



> C 



that in the anharmonic model we always have v 
at / — > 1, and although the simulation velocity is lower 
than v^ val ^ due to additional damping of the moving kink 
because of phonon radiation, the discrete kink still may 
reach a supersonic velocity. Thus, the variational ap- 
proach predicts that the supersonic kinks may be ex- 
pected in the anharmonic FK model only. 

Returning back to the discrete model, note that in the 
classical FK model the driven kink cannot reach even the 
sound velocity, because it exists some critical kink veloc- 
ity v c < c above which the driven kink becomes unstable 
and the system goes to the "running" state, where all 
atoms move with the velocity v ~ f jr\ [ ft6| . However, in 
the anharmonic FK model, the critical kink velocity may 
exceed the sound speed as has been observed already in 
the simulation |uj ]. The dependence of v c on the an- 
harmonicity parameter j3 is shown as inset in Fig. 2. In 
this calculation we used the following algorithm 111 ] : for 
a fixed value of / (we took / = 0.5), the friction was 
decreased starting from the overdamped case rj = 1 to 
the underdamped value n = 1CP 3 in 256 steps. At each 
step of rj decreasing we waited untill the steady state was 
reached and then checked if the transition to the running 
state took place. 



V. MULTIPLE KINKS 

To study multiple kinks with the help of a variational 
approach, let us consider the double kink as a sum of two 
single kinks separated by a distance r, 



u 2 (z) = u SG (z - r/2) + u SG (z + r/2). 



(18) 



Substituting the ansatz ( |18| ) into the functional (|11|), 
we obtain the effective energy E(d e g,r) which is now 
a function of two parameters d e g and r. Looking for 
a minimum of E(d c s,r) over d e s at r fixed, we found 
that for the classical FK model, a = 0, the function 
E(r) = min^jj E(d e g, r) is a monotonically decreasing 
function of r, i.e., the kinks are repelled from one an- 
other. On the other hand, for the anharmonic FK model 
the function E(r) has a minimum at some r = r min < oo, 
so two kinks attract one another and thus have to couple 
into a double kink. The "dissociation" energy of the dou- 
ble kink is very small at subsonic velocities, but becomes 
high enough at supersonic velocities (see inset in Fig. |^). 
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The "size" r m ; n of the double kink decreases with \v\ in- 
creasing. Thus, although the variational approach with 
the SG-type ansatz is too crude, it nevertheless predicts 
that the multiple kinks could be stable for high (super- 
sonic) kink velocities. 




FIG. 3. The same as in Fig. y for the double kink. Inset: 
the effective energy E(r) of the double kink as function of 
subkink's separation r for fixed kink velocities «k = 0, 0.9 c, 
and 1.8 c. 



both supersonic kinks and multiple kinks remain stable 
at nonzero system temperatures as well, at least for the 
time scale of our numerical simulation. 

Notice also that at high forces the kink velocity is close 
to that of the Toda soliton || . Indeed, the Toda soliton 
is characterized by the "jump" Au = 2[Ma//3, where ji is 
the parameter coupled with the soliton velocity v by the 
relationship v = csmh(jia) / [Mb. In the presence of the 
external substrate potential due to boundary conditions 
the jump Au must be equal to 2irp for the p-kink; thus 
we obtain fxa = 7r/3p, or 2/i = p(3. In particular, for the 
anharmonicity parameter (3 — 1/tt used in the simula- 
tion, we have v/c — (s'mh p)/p w 1.18 for the 27r-kink 
and v/c ~ 1.81 for the 47r-kink correspondingly. Thus, 
the supersonic and multiple kinks may be treated as Toda 
solitons "disturbed" by the external periodic potential. 



ACKNOWLEDGMENTS 

Discussions with A. R. Bishop, Yu. S. Kivshar, and M. 
Peyrard are gratefully acknowledged. This work was sup- 
ported in part by NATO Grant No. HTECH.LG 971372, 
INTAS Grant No. 97-31061, and by a grant from the 
Hong Kong Research Grants Council. 



As is well known, the SG kinks of the same topologi- 
cal charge always repel from one another. The same is 
true for static (v — 0) kinks of the discrete FK model, 
including the anharmonic (a ^ 0) model 0] (contrary 
to the variational approximation which mistakenly pre- 
dicted a weak attraction at \v\ < c). Thus, it must exist 
a threshold kink velocity v\ such that at small velocities 
< v < vi the steady-state solution of the model corre- 
sponds to the 27r-kink, while at high velocities v > vi, it 
corresponds to the double (4-7T-) kink (if v\ is lower than 
the kink velocity at / = 1, that is true at low enough 
values of rj). Indeed, the simulation results presented in 
Fig. |^ demonstrate that the double kink is stable within 
the force interval 0.5 < / < 0.8 but becomes unstable 
at higher as well as smaller forces, while the 27r-kink is 
stable for / < 0.5 only. Similarly one could expect the 
existence of a second threshold velocity v-i such that at 
v > t>2 the steady-state solution will correspond to the 
67r-kink, etc. 



VI. CONCLUSION 

Thus, we have shown that supersonic kinks as well as 
multiple kinks do exist in the driven discrete FK model, 
if the interatomic interaction is anharmonic. Note that 
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